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ABSTRACT 


An Isotropic sheet stiffened by means of an orthotropic quarter- 
plate Is considered. The quarter-plate Is assumed to be perfectly bonded 
to the metal sheet everywhere except an area of debonding which may de- 
velop due to high stress concentrations. The adhesive which has a small 
constant thickness will be treated as a shear spring. The loads are 
applied at Infinity and supposed to be transmitted through the metal 
plate. Shear stress distribution between the two plates are obtained 
from the continuity of displacements along an area where they are bonded 
to each ither. 


IN TRODUCTI ON 

As the composite materials are used In structures at an Increasing 
rate, the problems associated with the life prediction of such structural 
components become of major Interest. Among others, one can mention the 
use of such stiffening elements as stringers, strips and composite cov- 
er plates as some of the problem areas. Hence, significant number of 
solutions are now available In the open literature related to these areas. 

In this work, a sandwich panel will be considered where a metal 
sheet is assumed to be stiffened by a composite quarter-plate. Two 
plates are perfectly bonded to each other except an area of debonding 
which may develop during the life of the structure. Hence, the effect 
of the cover plate as a stiffening element and the way this Is influenced 
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by the debonding process will be Investigated. The loads are applied 
at the boundary of the Infinite metal sheet. 

The method of solution consists of satisfying the continuity of 
displacements In an area where two plates are perfectly bonded to each 
other [1], The solution of the resulting Integral equations gives the 
shear stress distribution between the two plates. Application of this 
approach requires the knowledge of Green's functions for both Isotropic 
([1]» [2]) and orthotropic ([3], [4]) plates which are listed In the 
Appendices. Although the solution Is presented for a metal sheet with- 
out a crack, by using the Green's function for a cracked plate given In 
Appendix A, results can also be obtained for a cracked sheet. As In- 
dicated In Appendix B, the quarter-plate solution Is based on an iter- 
ative process whereas the Green's function for an Infinite plate with 
or without a crack Is known In closed frrm ([1], [2] and Appendix A). 

Here, the boundary of the debond area will be assumed known. Other 
wise, it can be determined by an Iteration and using an appropriate 
failure criterion on the boundary [1], 

FORMULATION OF THE PROBLEM 

Let 0 denote the region where the two plates are bonded to each 

other and h c> h a and h p represent the thicknesses of the cover plate, 

the adhesive and the metal sheet respectively. The elastic constants 

for the orthotropic quarter-plate are E , E , G , and v and for the 

x y xy yx 

Isotropic plate are E and v. The shear modulus of the adhesive Is u . 
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T , T are the loads per unit length applied at Infinity. Due to the 
x y 

relatively small thicknesses of the plates and the adhesive, the problem 
will be solved under the assumption of generalized plane stress and the 
unknown shear stresses In the adhesive will be considered as body forces 
In the plate solutions [1]. This will satisfy the stress continuity 
conditions and the displacement conditions can then be written as follows: 

u (z) - u (z) » P(z) 
p c p a 

V (z) - v (z) - Q{z) , zcD . (1) 

where u . v and u , v are the displacements of the cover plate and the 
c c P P 

metal sheet In x,y directions, P and Q are th» shear stresses In the ad- 
hesive (see Fig. 2) and 

z * x + 1y (2) 


The displacements can be expressed as 


u p (x.y) 


1 


T x^y 

“V x + 


[k 1 (x,y,x 0 ,y 0 ) • 


P ( x o’ y o ) + k 2 (x ’ y * x o ,y o )Q( V y o^ dx o dy o + r ' b-d1sp1 


v p (x,y) = y + 4'nu MY+kT \ \ [k 3 (x,y * x o’ y o ) * 

K P P P n 


p ( x 0 - y 0 ) * ^ 4 (x,y,x 0 ,y 0 )Q(x 0 ,y 0 )] dx Q dy 0 + r-b-dlspl 


(3) 
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u c (*.y) s | | [k 5 (x *y* x o’y 0 )P(x o' y o ) + k 6 (x ‘ y * x o’ y o )Q(x o’ y o n 
D 

dx dy + r*b-d1spl . 
o o 

v c ( x ,y ) » | j [k 7 (x,y,x 0 ,y 0 )P(x 0 ,y 0 ) + k g (x,y,x 0 ,y 0 )Q(x 0 ,y 0 )] 

D 

dx dy„ + fb'dlspl . (4) 

o o 


where k«(3-v)/(1+v) and the kernels kj(x,y,x Q ,y 0 ) j s l,..4 are giver In 
Appendix A and kj(x,y,x 0 ,y 0 ) j=5...8 are obtained by an Iterative proc- 
ess In Appendix B [4].* 


Hence, from equations (l)-(4), the integral equations of the problem 
is obtained as 


P(x.y) + 


[[k n (x.y.x 0 ,y 0 )P(x 0 .y 0 )*k, 2 (x,y,x 0 ,y 0 )Q(x 0 ,y 0 )]dx 0 dy 0 


■ 'V'V EFT- x 

Pa 


Q(x.y) ♦ { { [k 21 (x,y,x 0 ,y 0 )P(x 0 ,y 0 ) + k 22 (x,y,x 0 ,y 0 )Q(x 0 ,y 0 )]dx 0 dy 0 
D 

s( VVi^^ • (x ‘ y)cD (5) 


Note that the location and orientation of the cover plate with respect 
to the metal sheet is insignificant if the metal sheet does not have a 
crack. In general, the rigid body displacements for both plates should 
be the same. 
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where 


k 11 “ ’ K7 ( Yk r k 5 ) 

a 

k 21 " “ h“ ^ k 3‘ k 7 ) 

d 

Y ■ 

The system of Fredholm integral equations (5) with logarithmic sing- 
ularities can be solved numerically as follows 

PUj.yj) ♦ y [ k ll (x j ,y j ,x ok ,y ok >Plx ok ,y ok ) * k l? (x j ,y J ,x ok ,y ok ) 

' Q(x ok- y ok )] Mk ' <V vT y> ETT^ x j 

Q(x j' y j’ * Jj [k 21 (x j ,y J ,x nk ,y ok )P( W y ok ) * k 22 (x j y J ,x ck ,y ok ) 

Q(x ok >y ok )1 Mk ■ ( V vT x> nnr y j " < 7 > 

p a 

where AAk's are the appropriately chosen area elements covering the en- 
tire domain D. Solving (7) the shear stresses P and Q at x 0(fl ,y 0k 
locations are obtained. 


k 12 " ' FT < yk 2- k 6> 

d 

P. 

k 22 * * IT < yk 4- x 8> 

d 

( 6 ) 
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APPENDIX A 


Kernels k J (x,y,x 0 ,y 0 ) J"l,..4: From [2] and [5] 

( a ) For a cracked Sheet [5] : 

k](x ( v.x 0 »y 0 ) « -Re [fj(z,z 0 ) ♦ f 2 (z.z 0 )] 
k 2 (x,y,x 0 ,y 0 ) • -Re [1^ (z,z Q ) - 1f 2 (z,z 0 )] 

k 3 (x y.x 0 .y 0 ) * -l n) [ f i( z ‘ z 0 ) + V z,z o^ 

^ 4 (x.y.x 0 ,y 0 ) - -J« [If^z^) - if ? (z ,z 0 )3 

(z ,Z 0 ) * -ic[log(z-z 0 ) + 1og(z-z 0 )] 

♦ | + [e 1 (!,j 0 ) + K 2 « 1 (i.i 0 

* (^1) [icfl 2 (z)-e 2 (z)]*(^2)e 6 (z,z 0 ) 

f 2 (z,z 0 ) * z~- - 1 + Ke 3 (z,z 0 ) - ° 3 (z» z 0 ) - 2 * z o^ + K " 4 ^ 2 ’ 


0,(z,zJ ■ log[z z -a 2 + /z*-a ? ] 


0 2 (z) ■ log(z + /z 2 -a 2 ) 


9 3 (z ’ z o ) * 


Z Q~ Z Q 

l/r^ 

0 


[1 + f(z,z Q )] 


(z-z) - . 

e.(z,zj = — f(z,Z ) 

4 0 2/j-JTa 7 0 


• 7 - 


(A.l) 


0 ) (A.2) 





■ jjgjf [«*•»•> - J <’o» 


(A. 3) 


I(z)-I(zJ 

* ° 


z-z. 


I ( z ) « /z ? -a r - z 


J(z) 


/Z*-fl2 


- 1 


z ■ x + iy 


z ■ x ? 1y„ 
o o •'o 


(A. 4) 


( b ) For a olaf* without a crack [2] : 


k 1 (x,y,x 0 ,y 0 ) - Z^e[log(z-z Q )] - Re[f 3 (z)] 


k 2 (x,y,x 0 ,y 0 ) ■ -Im[f 3 (z)] 
k 3 (x,y,x 0 ,y 0 ) » -Im[f 3 (z)] 


k 4 (x,y, x o,y 0 ) s &Re[log(z-z 0 )] + Re[f 3 (z)] 
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APPENDIX B 


Kernels kj(x,y,x 0 ,y 0 ) J-5...8: 

These kernels (Green's functions) can be obtained by an Iteration 
(successive approximation) scheme Introduced by Hetenyl and developed 
In [4] for orthotropic quarter-plate. 


Iteration starts with a basic symmetric system as Introduced In 
Fig. 3c. The stresses and the displacements for this system will be 
denoted by 


o5,oJ,Tjy,u c ,v c (x, y ) (B.l ) 

Then the solutions to two secondary half-plate problems are found 
(Figs. 3a and 3b) for R«1 as 


a a a a a, . 

°x’°y ,T xy ,u ,v ( x,y * 


and 


b b b b b/ « 
o x .ay,T X y.u .v (x,y) 


(B.2) 


which are Green's functions for half-plate problems. These can be ex- 
pressed either In r.»al variables [4] or In complex variables [5]. 


Due to the syntnetry, the basic system has zero shear stresses along 
the y axis, l.e., Tj y (0,y)=0, but oj^(0,y)»F o (y)^0. Hence, to obtain 
the solution to the quarter-plate, non-zero normal stress distribution 
F o (y) should be erased. Using the solution to the problem (3a) and 
applying a synmetric normal stress distribution o x (0,y)--F Q (y) on the 
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boundary we obtain the stress and displacement fields 

°x , V T xy ,u ' ,v (B.3) 

If this system Is superposed with the basic system, the resulting problem 
yields 


o x ( 0,y) « 0 but o y (x,0) ■ F^x) f 0 


(B.4) 


Hence second step will be to erase these non-zero o y (x,0) stresses a'ong 
the x axis. This can be done using the solution of problem (3b) thus 
yielding non-zero normal stresses along y axis again. Since In this 
approachjthe shear stresses along the x and y axes are always zero, 
repeating this process, one can obtain stress-free surfaces along these 
axes and the solution of the problem of a quarter-plate, within a reason- 
able degree of convergence [4]. The kernels (kg.k^) and kg, kg) can be 
obtained numerically by considering (P-1 • Q-0) and (P«0, Q»l) respect- 
ively. 


If one uses the expressions given In [4] for the solut 4 ons of prob- 
lems (3a) and (3b) the Iteration scheme will give 


o x (x,y) 


| k k ( k ) #co 

'(x,y) V - -1-^—1 — “f F 0 (n)dn[*i(x,y,n) ♦♦ 1 (x,y,-n)] (B.5) 

y v Vo 0 2 7 


T xy (x,y) 


2 

3 


where 
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(B.6) 


♦j(x.y.n) 

3 


<;bs 


\ /{Ckfx 2 My-n) 2 ][k|x 2 My-n) 2 ]) 


where 



?a 22 C2 a 1? ♦ *66 - /(2i 12* a 66 ,, ‘ 4# ll i 22 J 


a ll " b % * a 12 


W V.. . 1 

■-F-r • *»-r 

y x y 


a 66 " IT* 

xy 


(B.7) 


And the second iteration gives 


fo?(x,y) ] 


<J , 0 y(x.y) 

^yfx.y) 


k *k #°° 

- -A— F 1 (^)dC[^ 1 (x.y O^U.y.-O] 

1 r-0 $ 2 


(B.8) 


where 


(x-0 J y 
(x,y,C) -<jy* 

(x-dy 2 


[* /{[kj(x-C) 2 +y 2 ][kJ(x-C) 2+ y 1 ]} 


(B.9) 


etc. If these steps are repeated the final stress state can be obtained 
as 
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REJV.C 
OWCLnA i • 



° x («.y) 


o c 

X 

^° y (x,y) 

■ < 

°y \ 

T- tf (x,y) 
xy j 


k) 


(k,^k 2 ) 


t_k i k ? f CC^i (*.y .n) ♦ ♦i(x.y.-n)) 

1 ‘ Vo i \ 


«U .at 

l F m (n)]dn - [(^(x.y.U ♦ Vi(x.y.-C)) 

■"■0.2.4 C J 0 i { 


I. _ %(€)]«) 


m»l ,3,5 


(B. 10) 


And the functions F^.m-1,2.. are determined from the following 
recurrence relations 


iif+1 


Mo(ki*ko) f® 

• J - L = LJ - L '•“* 5F 


2xn a 

‘VMk’x’V) 


m'0,2,4. .. . 


F m f l (y) = " 2 L Fm(f)df 


ML 


c-o m (k}£y*)(lc^+y J ) 


m=l ,3,5. . . , ( B . 1 1 ) 

F 0 (n) Is obtained from the solution of problem (3c), [4]. 

This scheme Is applied to the displacements In an Identical manner 
by changing the functions $ and The use of complex variable solution 
[6] may prove to be much more convenient. 
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Fig. 1. Geometry of the problem. 





